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Abstract 

Coherent Harmonic Generation (CHG), and in particular Nonlinear Harmonic Gen- 
eration (NHG), is of importance for both short wavelength Free-Electron Lasers 
(FELs), in relation with the achievement of shorter wavelengths with a fixed 
electron-beam energy, and high-average power FEL resonators, in relation with 
destructive effects of higher harmonics radiation on mirrors, hi this paper we 
present a treatment of NHG from helical wigglers with particular emphasis on the 
second harmonic. Our study is based on an exact analytical solution of Maxwell's 
equations, derived with the help of a Green's function method. In particular, we 
demonstrate that nonlinear harmonic generation (NHG) from helical wigglers van- 
ishes on axis. Our conclusion is in open contrast with results in literature, that 
include a kinematical mistake in the description of the electron motion. 
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1 Introduction 



The study of Coherent Harmonic Generation (CHG) is of undisputed rel- 
evance in the field of Free-Electron Lasers (FELs). On the one hand, x-ray 
light sources based on Self- Amplified Spontaneous Emission (SASE) can 
benefit from CHG to radiate at shorter wavelengths at the same electron 
beam energy (see e.g. |jlt]). On the other hand, CHG can be a detrimental 
effect for high average power EEL oscillators, because of possible mirror 
damage from harmonics in the ultraviolet Thus, correct understanding 
of CHG is of interest for both short wavelengths and high average power 
applications, that are the two major up-to-date development paths in EEL 
physics. 

CHG is driven by bunching of the electron beam at harmonics of the fun- 
damental, and is characterized by the fact that harmonic components of 
the bunched beam radiate coherently [3]. The bunching mechanism may be 
linear or nonlinear. In the first case CHG is named Linear Harmonic Gen- 
eration (LHG). LHG arises when radiation at a certain harmonic induces 
electron beam bunching at the same harmonic. In the second case CHG 
is named Nonlinear Harmonic Generation (NHG). NHG arises when the 
intensity of the fundamental harmonicQ is strong enough to induce bunch- 
ing of electrons at different (higher) harmonics. In EEL processes, harmonic 
bunching of the electron beam due to interactions with the fundamental is 
always much stronger than that due to interactions with higher radiation 
harmonics. As a result only NHG has practical relevance, while LHG can 
be neglected. In general, CHG can be treated in terms of an electrodynam- 
ical problem where Maxwell's equations are solved with given sources in 
the space-frequency domain. Sources must still be obtained through the so- 
lution of self-consistent equations for electrons and fields. However, once 
these equations are solved one obtains macroscopic current and charge 
density distributions as a function of transverse and longitudinal coordi- 
nates at a given harmonic. Further on, solution of Maxwell's equations with 
these distributions as given sources characterizes harmonic radiation in the 
space-frequency domain. The dependence of sources in the space-frequency 
domain on transverse and longitudinal coordinates is complicated because 
is the result of the above-mentioned self-consistent process. However, here 
we deal with an EEL setup where an ultrarelativistic electron beam is sent, 
in free space, through an undulator with many periods. Then, paraxial and 
resonance approximation can be applied to simplify the characterization of 
CHG. In particular, for a fixed transverse position, the longitudinal depen- 
dence is always slow on the scale of an undulator period. 



Note that this argument is not restricted to the first harmonic. Here we have in 
mind an PEL lazing at the fundamental. 
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NHG has been dealt with in the case of a planar wiggler, both theoretically 
and experimentally, in a number of works Odd harmonics have 

maximal power on axi^ and are linearly polarized. Even harmonics have 
been shown to have vanishing on-axis power and to exhibit both horizontal 
and vertical polarization components. 

In this work we present the first exact theory of NHG from helical wig- 
glers. Our treatment is based on an exact solution of Maxwell's equations in 
the space-frequency domain based on a Green's function technique. In the 
next Section 2 electromagnetic sources are treated in all generality as given 
data to be obtained from self-consistent PEL codes. In particular, we will 
demonstrate that nonlinear harmonic generation (NHG) from helical wig- 
glers vanishes on axis. Later on we will focus on the second harmonic, and 
discuss a particular study-case. Our results are in contrast with conclusions 



in 11160, where NHG in a helical wiggler has also been addressed and the 
presence of on-axis power has been reported. In Section 3 we will show that 
the azimuthal resonance condition introduced in [16] to explain net energy 
exchange between particles and fields is a misconception. It arises from a 
kinematical mistake. Namely, the electron rotation angle in the helical wig- 
gler is confused with the azimuthal coordinate of the cylindrical reference 
system. Consequently, this misconception is passed on to simulations, lead- 
ing to incorrect results. We conclude our work with some final remarks in 
Section 4. 



2 Theory of nonlinear harmonic generation in helical wigglers 

2.1 Complete analysis of the harmonic generation mechanism 



In this Section we propose the first exact theory of NHG in helical wigglers. 
In particular we will give an analysis in the space-frequency domain of the 
second harmonic generation case. 

Paraxial Maxwell's equations in the space-frequency domain can be used 



to describe radiation from ultra-relativistic electrons (see |17l, uM). Let us 

call the transverse electric field in the space-frequency domain Ej_{z,rj_,co), 
where rj_ = xCx + yCy identifies a point on a transverse plane at longitudinal 
position z, Cx and Cy being unit vectors in the transverse x and y directions. 
Here the frequency o) is related to the wavelength A by cu = Inc/A, c being 
the speed of light in vacuum. From the paraxial approximation follows that 



^ Here we assume that the bunching wavefront is perpendicular to the (longitudi- 
nal) PEL axis. 
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the electric field envelope E± = E± exp [-icoz/c] does not vary much along z 
on the scale of the reduced wavelength A = A I {In). As a result, the following 
field equation holds: 



Ej_{z,fj_,co) 



(1) 



where the differential operator D is defined by 



2ico d \ 
c dzi 



(2) 



V^^ being the Laplacian operator over transverse cartesian coordinates. 
Eq. (1) is Maxwell's equation in paraxial approximation. The source-term 

vector f{z, r^_) is specified by the trajectory of the source electrons, and can be 

written in terms of the Fourier transform of the transverse current density, 

-± 

r_i, co), and of the charge density, p(z, r_i, co), as 



/= -471 



exp 



tcoz 
c 



(3) 



In this paper we will treat and p as macroscopic quantities, without 

-± 

investigating individual electron contributions. jj_ and p are regarded as 
given data, that can be obtained from any FEL code. Codes actually provide 
the charge density of the modulated electron beam in the time domain 
p(z, rx, f). A post-processor can then be used in order to perform the Fourier 
transform of p, that can always be presented as 



p = -p(z, fjL - rV(z)/ co) exp 



. So{z) 
ICO 



(4) 



where the minus sign on the right hand side is introduced for notational 
convenience only. Quantities r'o±{z), So(z) and Vg pertain a reference electron 
with nominal Lorentz factor yo that is injected on axis with no deflection 
and is guided by the helical undulator field. Such electron follows a helical 
trajectory r'o±{z) = r'^^Cx + '^'oy^y We assume that the reference electron rotates 
anti-clockwise in the judgement of an observer located after the undulator 
and looking towards the device, so that 



Kxiz) = —r[(^os{ku,z) - 1] , r'{z) = 



K 



sin(fca,z) . 



(5) 
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In Eq. (5), K = AivCH^/ (InmeC^) is the undulator parameter, A^u = 2n/kzu being 
the undulator period, (-e) the negative electron charge, H^y the maximal 
modulus of the undulator magnetic field on-axis, and the rest mass of 
the electron. The corresponding velocity is described by Vo±iz) = VoxSx + Vgy^y 
with 



Kc Kc 
Vox{z) = sin(fc,„z) , Voyiz) = — cos(fc^„z) . (6) 

Yo To 



Finally, So(z) is the curvilinear abscissa measured along the trajectory of the 
reference particle. 

One may always present p as in Eq. (4). However, introduction of p is 
useful when p is a slowly varying function of z on the wavelength scale. 
This property is granted by the fact that the charge density distribution 
under study originates from an EEL process. Erom this fact it also follows 
that p is slowly varying on the scale of the undulator period Aw and is 
peaked around each harmonic of the fundamental cOr = Ik^cfl, that is fixed 
imposing resonance condition between electric field and reference particle. 
The word "peaked" means that the bandwidth of each harmonic component 
obeys Aco/Qiojr) « 1 for each positive integer value /z. Here = l/il-vl^/c^) 
is the longitudinal Lorentz factor. Note that, for the reference electron, fz 
does not depend on z. We have 



1 y 



where the last equality follows from v^^ = ~ ^o±' together with Eq. (6). 
Einally, the relative deviation of the particles energy from yotneC^ is small, 
i.e. 6y/yo «: 1. 

In this paper we will be interested in the case when the transverse beam 
dimension Oj_ is much larger than the electron rotation radius r^i,, i.e. Oj_ » r^- 
We might then substitute the dependence on rj_ - r'o±iz) in Eq. (4) with a 
dependence on rj_, because the electron rotation radius is negligible with 
respect to Oj_ and individual electrons can be considered as occupying a 
fixed transverse position. However, we will not do so. In fact, starting from 
Eq. (4) to develop our theory, we will be able to crosscheck our results for an 
extended source with the well-known asymptotic for a filament beara (see 



Section 2.2), i.e. Oj_ « r-^, and demonstrate agreement with results in |20|]. 
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We note that for a generic motion one has 



^s(Z2)-s(Zi) Z2-Z1 



CO 



2yl{z)c 



(8) 



Thus, with the help of Eq. (4), Eq. (3) can be presented aS 



f = 4n exp 



z 

I 



dz- 



CO 



2ylc 



ptz;fl -rV(z)] , 



(9) 



where we used the fact that = Vo±p- In fact, for each particle in the 

beam 6y/yo ^ 1. Therefore we can neglect differences between the average 

-± 

transverse velocity of electrons {Vj_) and Vo±, so that ij_ = {Vj_)p - Vo±P- 

We will now introduce a coherent deflection angle if'^\ where the superscript 
(c) stands for "coherent", to describe the transverse deflection of the electron 
beam as a wholS This means that we account for a possible deflection 
angle ff'^^ in the the trajectory of reference electron. We therefore perform 
the following substitutions in Eq. (9): 

rV(z)^rf (z,rf^') = 

rV(z) + ff'h = 1^ [cos(fc,„z) - 1] + r]i'^z| + |^ sin(fc^z) + i]fz^ Cy , 

(10) 



Vo±{z) + af'^ 



Kc . . 
sm{h,z) + CI] 

y 



— cos(/Cj„z) + cr]\. 

y ^ 



(c) 



(11) 



Using v±{z, rf'^'>) in place of Voi.{z) implies that yz{z, rf'^^) is now a function of 
both z and rf'^\ In particular, l/yz(z, rf '^^) = 1 - vl{z, rf'^^)/c'^, where vl = v'^- v]_ 



Eq. (9) is the analogous of the source term on the right hand side of Eq. (11) in 
reference [15], dealing with planar undulators and an electron beam modulated at 
the second harmonic only. There is only a slight notational difference in that the 
symbol jO in the present paper corresponds to jgai/c in reference listl . 
^ With this, we assume that the deflection angle ff'^^ is constant. This is the case only 
if we do not account for focusing elements within the undulator. Generalization to 
account for betatron motion of electrons is, however, straightforward. 
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is the square of the electron longitudinal velocity. It follows that l/ylin Eq. 
(9) should also be substituted according to 



1 



(^^1' + {4^f] + y [-^x ' sin(fc,„z) + r^f cosfcz)] , 



(12) 



yielding 



f = 4n exp 



i dz 



CO 



lcyl{z,Tf^)) 



'^Mz, f^) - V, p[z- - (z, rf^^)\ 



(13) 



We find an exact solution of Eq. (2) without any other assumption about 
the parameters of the problem. A Green's function for Eq. (2), namely the 
solution corresponding to the unit point source can be written as (see |ll5|): 



G(z„-z';r--rl) = -^^^^^exp 

assuming Zo - z' > 0. When Zo - z' < the paraxial approximation does not 
hold, and the paraxial wave equation Eq. (1) should be substituted, in the 
space-frequency domain, by a more general Helmholtz equation. However, 
the radiation formation length for Zo - z' < is very short with respect to 
the case Zg-z' > 0, i.e. there is no radiation for observer positions Zo - z' < 0. 
As a result, in this paper we will consider only Zg-z' > 0. It follows that the 
observer is located downstream of the sources. 

This leads to the solution 



ICO- 



'1.0 



2c{Zo - z') 



(14) 



EAzo,Ao) = -J dz' ^ ^ df 

— CO 

xp{z' ,r'^ - j*;^(z',7](*^))) exp 



'■^vdz',Tf^^)-y'^ 



ICO 



2ciZo - z') 



+ i I dz 



CO 



'2cy2(z,rf^)) ' 



(15) 



—f 

where represents the gradient operator (acting on transverse coordinates) 
with respect to the source point, while (Zo,rj_o) indicates the observation 
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point. Note that the integration is taken from -oo to oo because the sources 
are understood to be localized in a finite region of space (i.e. p is different 
from zero in a finite region of space). This is in agreement with Zg-z' > and 
within the applicability criteria of the paraxial approximation. Integration 
by parts of the gradient term leads to 



oo 

C J Zo-Z' J 



C Zq-Z' 

Xp (z', r'x. - f^l\z', exp [/Ot(z'/ r'±, r/W)] , 



(16) 



where the total phase Or is given by 
2c{Za-z')\ J 



(^T — CO 



dz 



CO 



2c72(z,rf^)) 



(17) 



We now make use of a new integration variable /' = r'± - f[\z', J]^''^) so that 



00 

t = J^ fdz'^ [df 

C J Zo-Z' J 



(vAz'.nh no-^^{z\iih-i 



n 



Zo-Z' 



xp(z'J)exp[i^T{z',lr]^] , 



(18) 



and 



^o^^?(zVr^-/'l 
2c(Zo - z') 



/|2 



+ 



z 

/ 



dz 



Ci) 



lcy\{z,vt'~^) 



(19) 



We will be interested in the total power emitted and in the directivity dia- 
gram of the radiation in the far zone. 

We therefore introduce the far zone approximation calling the observation 
angle Q = r^olzo, setting 6 = \d\ and taking the limit for Zo » L^, where 
= ^w^w is the undulator length: 
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X exp 



z 

^ [zod^ - 2d- (z',T]W) -26-1'+ z:&\ + i j 



dz 



CO 



2cyl{z,rf^)) 



(20) 



Substitution of Eq. (10), Eq. (11) and Eq. (12) in Eq. (20) yields the following 
field contribution calculated along the undulator: 



Ex — — I dV exp [zOo] ^ dz'p[z',l') exp [zOr] 

X 1^ sin {Kz') + (d, - rji^') j + cos (Kz') + (dy - r]ffj 



(?1) 



where 



O7 



coz 
'2c 



Key 

cyky, 



z 

[{Oy - r]f) sinih^z') + (0, - 1]'^^) cos{h,z')] 



(22) 



and 



2 Ky 



- (0,/; + OyVy) 



(23) 



Note that the integration in Eq. (21) is performed in dz' over the undulator 
length, i.e. is limited to the interval [-Lwl2, Lf„/2]. The reason for this is that, 
working under the resonance approximation in the limit for Ny, » 1, one 
can neglect contributions to the field due to non-resonant elements outside 
the undulator 0]. 

We are interested in studying frequency near the fundamental harmonic 
cOr = 2ki„cyl or its h-th integer multiple. We specify "how near" the frequency 
CO is to the h-th harmonic by defining a detuning parameter C;, as 



Ch = 



CO 



2yic 



hk -—k 

COr 



(24) 



Here 00 = hcor + Acu. Eq. (22) can thus be rewritten as 
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Key 



so that altogether, Eq. (21) can be presented as: 



(25) 



£, = ^ J rf/; J di; J d2:p[2:,vj^ 

-oo -oo -L„,/2 



X exp ^ I— 



X exp [i \hk^ + Q + g (0, - r]f )' + 1^ (0y - i^O'] 2' 



/-^ [(0, - 1]^^) siniKz') + (0, - 4^)) cos(?c,„z')] 



X 



+ 



^ (exp[zfc^„z'] - exp[-zfc^„z']) + (0;, - 7]if') 

~2y (^^P[*«'^'] + exp[-zlj„z']) + (0^ - r]f) 
We make use of the well-known expansion (see [19]) 

oo 

exp [ia sin (xp)] = ^ /p(fl) exp [ipip] , 



(26) 



(27) 



where }p indicates the Bessel function of the first kind of order p, while a 
and \p are real numbers. Eq. (26) may thus be written as 



00 oo 



Lu,/2 



'^Jdl'.j dl'y j dz'p[z',l'J'y) 

-Lu,/2 

Zo(0? + 0;) m - rii^)) 



-oo —oo 



.CO 

X exp ^ I— 



X ^'^P ^ [^('^ + + ^)fczi;Z'] 

m,n=-oo 



X 



^ (exp[zlj„z'] - exp[-zfc^„z']) + (0;, - rf^'^) 
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+ 



— (exp[iku,z'] + exp[-zla,z']) + [dy - rfy^) 



27 



where 



u = ; and V = — 



cyku 



cyku 



(28) 



(29) 



Whenever 



CO 



(30) 



the second phase factor in z' in Eq. (28) (the one containing C/,) is varying 
slowly on the scale of the undulator period A-a,. As a result, simplifications 
arise when N^, » 1, because fast oscillating terms in powers of exp[ikji,z'] 
effectively average to zero. 

Let us use Eq. (28), to discuss harmonic radiation characteristics. Within the 
resonance approximation we further select frequencies such that 



|Aa;| 



COr 



«: 1 , i.e. IQI <c k^ . 



(31) 



Note that this condition on frequencies automatically selects observation 

angles of interest h{d - if'^^f- «; l/y^. In fact, if one considers observation 

angles outside this range, condition (30) is not fulfilled, and the integrand in 

Eq. (28) exhibits fast oscillations on the integration scale L^. As a result, one 

-* 

obtains zero transverse field, = 0, with accuracy 1 /Ni„. Under the con- 
straint imposed by (31), independently of the value of K and for observation 
angles of interest we have 



(32) 



and similarly 



1 + ^\ yMOy - r//^'| < yMOy - « 1 • 



(33) 



This means that, independently of K, \v\ <c 1 and \u\ «; 1, so that we may 
expand ]n{v) and /^(m) in Eq. (28) according to /p(x) ^ SP(p)[2-IPl/r(l + \p\)] jIpI, 
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where S(p) gives the sign of the integer pifp ^0 and unity if p = 0, while T 
is the Euler gamma function 



r(z) = j ^'^ exp[-f] . (34) 



From now on we will restrict our attention to the second harmonic, noting 
that we may study any harmonic of interest in a similar way. 

Thus, h = l. Terms giving a non-zero contribution after integration in dz' in 
Eq. (28) are those for n = -m -3,n = -m - 1 and n = -m - 2. Therefore, we 
rewrite Eq. (28) as: 



oo oo Lii,/2 f 



-Ln,/2 



xp (z', /;) exp {i [C. + 1^ (0. - + 1^ (0, - r/i^)f ] z'} 



X 



+ 



exp 



imn 



-(/- 
2i7 V 



2 J 

m=— oo 

(Qx - r]x^) }-m-2{v) exp [-in] 



m-3 



(v) exp 



i3n 
~2~J 



/-m-l(u)exp -y j 



gr + 



m-3 



(u) exp 



3m 
'~2~J 



-/-„-i(t;)exp|-y 



m 



j + (0y - ^ly') /-m-2(^^) exp [-Z7l] 



(35) 



where we neglected the phase contribution K{dx- rfx)IQ^wy) because is 
much smaller than unity. Non-negligible terms in the expansion of /p(-) 
are those for small values of \p\, because arguments are much smaller than 
unity. 

Let us begin with contributions for m = 0. In this case /^(m) = /o(m) ~ 1 and 
}-fn-i{v) = }-i{v) ~ V. Thus, terms in }-m-i{v) scale as Kv/y. Terms in /_m-3(u) 
scale as Kv^/y and are negligible. Terms in J-m-2{v) scale as \dx,y - i{x!y\ ~ 

^ 2 

^Kvy^ - rf'^) / Vl + «; ¥\v\ly and are also negligible. 

In the case m = -1, fm{u) = J-i{u) ~ u and f-m-i{v) = Jo{v) ~ 1. Thus, terms in 
f-m-i (^) scale as Ku / y, that is of the same order of Kv jy. Terms in /-m-s (^) scale 
as Kuv^/y and are negligible. Terms in /_m_2(t') scale as \6x,y - rj^^yl \u\ \v\ ~ 

2 

4Kuy^ d-ff'^ / Vl + X2 « K\u\/y and are also negligible. 

Similarly, all other values of m give negligible contributions. As a result we 
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obtain 



(^x + ^j) 



X J dl'^ J dl'y J dz'expli^ 

-OO -OO ^ 

xp (z', /;) exp [iC2z'] exp {f |^ [(0, - r^'^^f + (Oy - r^f) 



(36) 



Note that the electric field is left circularly polarized (rotating anti-clockwise 
in the judgement of an observer located after the undulator and looking 

towards the device, as it must be) and vanishes at = rf'^\ Polarization 
characteristics are the same as for the fundamental harmonic, although the 
fundamental does not vanish at = rf-'^'. Spatial resonance is organized 
along the undulator for particular values of m, as discussed above. 

We conclude our analysis of NHG in helical wigglers studying on-axis 
harmonic generation. We can do so in all generality, i.e. for any harmonic 
number, with the help of Eq. (28). We set O-rf^^ = 0, because the concept of 
"on-axis emission" is evidently related to the coherent deflection angle rf'^\ 
We see by inspection that when - ff'^^ = 0, the only non-zero contributions 
are for n = w = 0, because both u - and v = 0. Thus, Eq. (28) can be 
rewritten as 



Lw/2 



E^ = !^ J df J dz'p(z'j)exp[iCHz'] 



-U,I2 



X ] I -^(exp[z(/z + l)Kz'] - exp[z(/z - l)Kz"[) 



+ 



liy 
K 



— (exp[z(/z + l)Kz'] + exp[i{h - l)k^z']) 



27 



(37) 



Since jo is a slowly function of z' on the scale of the undulator period and 
Ch ^ K>, we see by inspection that, after integration in dz', one obtains 
non-zero on-axis field only for /z = 1, leading to 



£x = - 



Ka){x + iCy) 
Iczg-y 



hJ2 



Jdf J dz'p{z'J')exp[iChz'] . 



(38) 



All other harmonics vanish on-axis. This result is in open contrast with what 
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reported in reference [|16|1 . In Section 3 we will elaborate further on this issue. 



2.2 Analysis of a simple model 



Let us treat a particular case to exemplify our results. Namely, let us consider 
the case when C2 = (i.e. co = IcOr) and p {z' , l'^, l'^ is given by 



( 7'2 



P 



Inco 



exp 



^ y 



2\ 



exp 



20L>r 



[4X + rifly) 



HlJz) , 



(39) 



with Hl^„ a window function equal unity inside the undulator and zero 
everywhere else. Here Ig is the bunch current, is a constant determining 
the strength of the bunching and cr^ the rms transverse size of the electron 
beam. 

This particular case corresponds to a modulation wavefront perpendicular 
to the direction of motion of the beam|£l In this case Eq. (36) can be written 
as 



E, = 



ioj^lo^i {e^c + iey) 

nolc^Zo l + K^ 



[(e.-r^^yi(e,-nf)]exp 



L c 



i-Zo{el + el) 



00 oo 



X 



j dv, j di[j j dz' exp [(0, - 4^) /; + [Oy - nf) i;]}h,jz') 

—00 

(e,-,,?'f + (e,-,,<f 



.COr 

X exp \ I— 



z fexp 



^ y 



(40) 



Direct calculations yield 



Information about the modulation wavefront of the beam is included in p. A 



phase factor of the form exp 



describes a plane wavefront. The case 



icoa T]'"^) ■ r /c 

a = 1 that is considered here corresponds to a modulation wavefront orthogonal to 
the z axis. When a = the modulation wavefront is orthogonal to the direction of 
propagation. Other values correspond to a modulation wavefront that is orthogonal 
neither to the z axis nor to the direction of propagation. This may be the case in SASE 
XFEL setups with very long saturation lengths (order of lO'^ m), where there can be 
coherent orbit perturbations. Note that if the beam is prepared in a different way 
so that, for instance, the modulation wavefront is not orthogonal to the direction of 
propagation of the beam, Eq. (36) retains its validity. 
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X exp 



xsinc 
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( 2o^, col 
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2 



[(e,-^?) + i(e, -,<")] 
(e,-4f + (e,-'(?y 



(e,-,?)%(e,-,l,'')' 



(41) 



Going back to our particular case in Eq. (41), a subject of particular interest 
is the angular distribution of the radiation intensity, which will be denoted 
with I2. Upon introduction of normalized quantities 



d 



I IcOrL-a, 



x,y ■ 



ilx,y 



and of the Fresnel number 



N 



cLu 



(42) 



(43) 



one obtains 



(|l-f|) 



I = const X 



6 -f] exp < -N 6 - i] 



^c) 



smc - 



0-77 



(44) 



In the limit for N ^ 1, Eq. (44) gives back the directivity diagram for the 
second harmonic radiation from a single particle in agreement with [ jiol ]. 

As an example, the directivity diagram in Eq. (44) is plotted in Fig. 1 for 



different values of N as a function of 
intensity J^^" at each value of N. 



e-f\ 



(c) 



, normalized to the maximal 



The next step is the calculation of the second harmonic power that is given 
by 



W2 = ^J" dXo J dyo\E^{Zo,Xo,yo,tW = ^ ^ dXo J dyo\E^{Zo,Xo,yo)? , 

— 00 —CXD —00 —00 

(45) 
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Fig. 1. Plot of the directivity diagram for the radiation intensity as a function of 

:? -»(c) 

6 -T] for different values of N, normalized to the maximal intensity at each 
value of N. 



where (...) denotes averaging over a cycle of oscillation of the carrier wave 

-> 

and Ej_{Zo,Xo, yo, t) is the electric field in the space-time domain at position 
{Xo, yo, Zo) and time f . 

We will still consider the model specified by Eq. (39) with C2 = 0. It is 
convenient to present the expressions for W2 in a dimensionless form. After 
appropriate normalization it is a function of one dimensionless parameter 
only, that is 



= F2(N) = In (1 + ^). (46) 



Here W2 = ^Nil^l is the normalized power, while the normalization con- 
stant is given by 



(47) 




For practical purposes it is convenient to express Eq. (47) in the form 

^^f' = {^^ ^^^4^] ' (48) 



^1 + ^2/ ^Xyh 
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Fig. 2. Illustration of the behavior of F2{N). 

where Wi, = mgC^ylg/e is the total power of the electron beam and I a = 
meC" le - 17 kA is the Alfven current. 

The function FiiN) is plotted in Fig. 2. The logarithmic divergence in F2(N) in 
the limit for N «: 1 imposes a limit on the meaningful values of N. On the one 
hand, the characteristic angle ©max associated with the intensity distribution 
is given by 0^^^ ~ 1/N. On the other hand, the expansion of the Bessel 
function performed in Eq. (36) is valid only as 6^ < Nj„. As a result we find 
that Eq. (46) is valid only up to values of N such that N > N~^. However, in 
the case N < we deal with a situation where the dimensionless problem 
parameter N is smaller than the accuracy of the resonance approximation 
~ N~^. In this situation our electrodynamic description does not distinguish 
anymore between a beam with finite transverse size and a point-like particle 
and, for estimations, we should replace In (N) with In {N~^). 



3 Discussion 

NHG in a helical wiggler has been addressed in In the numerical 
study-case proposed in that reference, an ultrarelativistic beam at 140 MeV 
is considered, driving an EEL oscillator operating near 1 /^m wavelength 
in free space. NHG is studied through simulations for an undulator with a 
uniform field region of 20 periods. The complete parameter set can be found 
in reference [16]. Eor our purpose, parameters reported above are enough 
to guarantee that paraxial and resonance approximation can be applied to 
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the study-case in |16|1 . In that reference one can read: "Our conjecture was 
that on-axis harmonic excitation due to NHG in helical wigglers can arise 
because the nonlinear bunching due to the fundamental creates on-axis 
harmonic radiation. This conjecture is borne out in simulation", and also: 
"Because NHG is driven by the fundamental, which excites on-axis modes, 
we speculate that NHG in helical wigglers will have substantial on-axis 
power". In the same reference one finds a statement regarding harmonic 
radiation from a single electron too: "Kincaid |20J showed that spontaneous 
generation from helical wigglers vanishes on axis". 

Thus, the main result of (l^ is that characteristics of helical undulator radi- 
ation from an extended source, i.e. a bunched electron beam, are drastically 
different compared to those from a single electron studied in [20]. In the 
previous Section we have seen that our conclusions are in open contrast to 



1 16]. In this Section we give reasons why, in our understanding, results in 
|16|] are incorrect. 



A first hint follows from general properties of linear superposition. Any lin- 
ear superposition of a given field harmonic from single electrons conserves 
single-particle characteristics like parametric dependence on undulator pa- 
rameters and polarization. In particular, if a field harmonic from a single 
electron vanishes on-axis, it must vanish for the linear superposition as well. 
Consider a bunched beam prepared in such a way that electrons enter the 
undulator with fixed offsets with respect to the longitudinal axis and with 
specific phases related with their positions along the bunch. Radiation fields 
generated by this beam can be seen as a linear superposition of fields from 
individual electrons. Now, as shown in [20], harmonic radiation from each 
of these electrons vanishes on axis. It follows that the total field has zero 
on-axis power as well. This property is conserved when the dependence of 
charge and current density distributions of the sources on the longitudinal 
coordinate is slow on the scale of the undulator period. This is always the 
case for NHG from ultrarelativistic beams in an PEL system in free space. 



This argument suggests that results in reference [|16l] are incorrect. A formal 



demonstration of this fact can be given discussing the mechanism proposed 



in 11160 to explain net energy exchange between individual electrons and har- 
monics of the fundamental in helical undulators. This is obtained through 
an azimuthal resonance condition, that is introduced with these wordilf 
"The azimuthal electron motion in helical wigglers is = k^z {kw is the 
wave number for the wiggler period A^o), which couples to circularly polar- 
ized waves that vary as exp(z(^;,), where (/);, = kz + hd - cot is the wave phase. 



^ Note that this mechanism is independent of the type of harmonic generation 
process, because it involves interaction between a single particle and the electro- 
magnetic field. In particular, it pertains both LHG and NHG. 
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Fig. 3. Cylindrical coordinate system and undulator setup. 

Hence, the phase along the particle trajectories varies as (ph = (k + hk^o)z - cot, 
and the h\h. order azimuthal mode corresponds to the h\h. harmonic reso- 
nance [i.e., cu ~ (/c + hk^)Vz]" (cited from Here indicates an azimuthal 
position as in Fig. 3. 

It should be remarked that for h = 1 one obtains co x (k + ku)Vz that is the 
usual resonance condition between an azimuthal symmetric wave and an 
electron in a helical undulator. This means that the phase (/);, = kz + hd - cot 
pertains a circularly polarized wave whose electric field is written in terms 
of unit vectors Bp and ee in polar coordinates, and not in terms of unit vectors 
Cx and Cy in cartesian coordinates (Cp, Cq, and Cy are defined in Fig. 3). In 
fact one may write the electric field of a (e.g. left) circularly polarized plane 
wave at position r and time t as : 



£(r, = Eo{ex + i Cy) exp 
= Eo{er + i ee) exp 



ikz + i{h - 1)6 - icot^ 

ikz + ihd - icot^ = Eo{er + i e^^ exp ]^(ph^ , (49) 



Eo being a constant field strength. This definition of (ph is in agreement with 
notation in [|21, 22], and justify words in |16]: "the h\h. order azimuthal mode 
corresponds to the /ith harmonic resonance". In fact, we can write the angle 
\\) between the transverse velocity of the electron and the transverse electric 
field vector as i/^ = kj„z + [kz + {h- 1)6 - cot]. Authors of |ll6i|2l|] use relation 
6 ~ k^aZ in the expression for cph, and automatically in the expression for 
\p too. Using also dz = Vzdt, where Vz — c is constant for an electron in a 
helical undulator, they obtain resonance for d\p/{dz) = k + hk^a -co/Vz = 0. As 
said before, for the particular case h = 1 (fundamental harmonic) one has 
h - 1 = (azimuthal-symmetric wave) and dip/{dz) = for co = (k + k-u,)Vz, 
that is the usual resonance condition. 



We argue that the azimuthal resonance described above is a misconception 
arising from a kinematical mistake. Namely, it is incorrect to use relation 
6 » k^„z in the expression for xp as done in [|16l,l21l] . Let us show this fact. 

We begin introducing some characteristic scale of interest. On the one hand. 
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electron rotation radius ~ r„ electron beam area ~ cr| 
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^ / 
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single electron diffraction area ~ XL^ 

6 = constant 



Fig. 4. Hierarchy of characteristic scales of interest. 



the radiation diffraction size for a single particle is of order V^L^. On the 
other hand, the electron rotation radius is given by Tyj = {K/y)Azi,. It follows 
that 



' w 

AL,, 



^^A^ = 



1 2K^ 
AU„ 1 + 



AAu, = 



«1, 



(50) 



where we used the fact that 7/^ = (1 + K^)yl, where fz has been introduced 
in Eq. (7). Inequality (50) holds independently of the value of K, because 
Nf„ » 1. Thus, the electron rotation radius is always much smaller than 
the radiation diffraction size. To complete the picture we introduce a last 
characteristic scale of interest in our problem, pertaining the electron beam 
rather than a single electron. This scale is the transverse size of the electron 
beam, Oj_. Straightforward geometrical considerations show that it makes 
sense to talk about aj_ only when o]_ » r^^. This is the case in practical 
situations of interest. There is still room to compare the beam size Oj_ with 
the radiation diffraction size. In the case Oj_ <^c V^X^ we deal with a filament 
electron beam, and results in [20] can be directly applied. In the opposite 
case the filament beam approximation breaks down. We are interested in 
this last case, where single-particle results cannot be used. We have therefore 
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established a hierarchy in the characteristic scales of interest: cr^ > AL^y » 
The situation is described in Fig. 4, where we schematically indicated the 
electron beam area and the diffraction area as disks. In Fig. 4 we also indicate 
the azimuthal and radial coordinates of a single electron, respectively 6 
and r. As one can see, the azimuthal coordinate of each electron, 6, is fixed 
during the motion inside the undulator with the accuracy of the resonance 
approximation, scaling as l/N-a,. In contrast to this, the identification 6 » kiuZ 
is made in [|l6l, l2l|l . This is a kinematical mistake. If ~ fcj^z each electron 
would be rotating around the origin of the coordinate system, that is not 
the case. Thus, we have shown that ip = k^„z + [kz + {h - 1)6 - cot] with 
6 constant for each electron along the trajectory in the undulator. As a 
result d\p/{dz) = only for h = 1 that yields the usual resonance condition 
CO = {k + k^)Vz. Summing up, our conclusion is that a kinematical mistake is 
at the origin of the azimuthal resonance described in |16l,l2l|l . The azimuthal 
resonance condition is a misconception following from this mistake. This 
misconception is subsequently passed on to simulations in |16|l, resulting 
in incorrect outcomes. Harmonic emission exists for a single electron [20], 
and it also exists for an electron bunch. However, qualitative properties are 
different with respect to what has been predicted in [|l6l] . In particular, as 
we have seen in Section 2, on-axis power vanishes. 



4 Conclusions 



In this paper we discussed NHG in helical undulators, with particular em- 
phasis on second harmonic generation. First we discussed the NHG mech- 
anism in helical undulators in all generality. Then we specialized our study 
to the case of second harmonic generation. Finally, to exemplify our results, 
we treated a simplified model where the beam modulation wavefront is or- 
thogonal to the z axis, has a Gaussian transverse profile and is independent 
on the position inside the undulator. 

Our results show that on-axis harmonic generation from helical wigglers 
vanishes. This applies to any harmonic of interest with the exclusion of the 
fundamental, and independently of the form of the electron beam mod- 
ulation (assuming that the electron beam as a whole propagates on-axis). 
Important consequences follow regarding the two mainstream development 
paths in FEL physics. First, as concerns short wavelength (x-ray) SASE EEL 
devices, vanishing on-axis harmonics make the option of a helical undula- 
tor less attractive as regards the exploitation of NHG radiation. Second, as 
concerns high average-power FEL oscillators, vanishing on-axis harmon- 
ics suggest that helical undulators carry relevant advantages over planar 
undulators, as potential for mirror damage is reduced. 
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Previous studies reported non-vanishing on-axis power, due to the nature 
of a particular azimuthal resonance condition. We showed that this res- 
onance condition is a misconception, arising from a kinematical mistake. 
This misconception was passed on to simulations, that confirmed the pres- 
ence of on-axis power out of NHG from helical wigglers. This result is also 
incorrect. Simulations are undoubtedly very important scientific tools, but 
they follow specific models. The correctness of their outcomes is related to 
the correctness of these models, meaning that validity of simulations should 
always be cross-checked with analytical results from simplified study-cases. 
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